Optimal Choices of Reference for Quasi-local Energy 
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We have proposed a program for determining the reference for the quasi-local energy defined in 
the covariant Hamiltonian formalism. Our program has been tested by applying it to the spherically 
symmetric spacetimes. With respect to different observers we found that the quasi-local energy can 
be positive, zero, or even negative. The observer measuring the maximum energy was identified; the 
associated energy is positive for both the Schwarzschild and the Friedmann-Lemaftre-Robertson- 
Walker spacetimes. 
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Introduction. An outstanding fundamental problem in 
general relativity is that there is no proper definition 
for the energy density of gravitating systems. (This 
can be understood as a consequence of the equivalence 
principle.) The modern concept is that gravitational 
energy should be non-local, more precisely quasi-local, 
i.e., it should be associated with a closed two-surface 
(for a comprehensive review see Here we consider 

one proposal based on the covariant Hamiltonian formal- 
ism [2J wherein the quasi-local energy is determined by 
the Hamiltonian boundary term. For a specific space- 
time displacement vector field on the boundary of a re- 
gion (which can be associated with the observer), the 
associated quasi-local energy depends not only on the 
dynamical values of the fields on the boundary but also 
on the choice of reference values for these fields. Thus a 
principal issue in this formalism is the proper choice of 
reference spacetime for a given observer. 

It is generally accepted that the gravitational energy 
for an asymptotically flat system should be non-negative 
and should vanish only for Minkowski space (see, e.g., 
for proofs of this property for GR see [1, Hj]). In view of 
this it has been natural to regard these properties as de- 
sirable for a good quasi-local energy [l|, Q • The idea that 
a suitable reference should be the one which gives the 
minimal energy then followed quite naturally. We have 
proposed using this approach to choose the optimal refer- 
ence for the covariant Hamiltonian boundary term. Here 
we consider this optimal choice program for both static 
and dynamic spherically symmetric spacetimes (the most 
important test cases). Specifically, for the Schwarzschild 
and the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
spacetimes, we found the resultant quasi-local energies 
can be positive, zero, or even negative for different ob- 
servers. However, for both cases, there is one observer 
who would measure the maximum energy, and for this 
observer the associated energy is positive. Furthermore 
we find that this energy-extremization program (at least 
for these spherically symmetric systems) is equivalent 



to matching the geometry at the two-sphere boundary, 
which provides for a simple interpretation of the displace- 
ment vector. 

The Hamiltonian Formulation. We begin with a brief re- 
view of the covariant Hamiltonian formalism 0, H, @, [13] • 
A first order Lagrangian 4-form for a /c-form field ip can 
be expressed as C — dip hp — A(ip,p). The action should 
be invariant under local diffeomorphisms, which infinites- 
imally correspond to a displacement along some vector 
field N. From Noether's theorem there is a conserved 
translational current which can be written as a 3-form 
linear in the displacement vector plus a total differen- 
tial: H(N) := £ N ipAp-i N £ := N^H^ + dB(N). Here 
H„ = -i„ip A f| + A i M p with ? := (-l) fe ; this iden- 
tity is a necessary consequence of local diffeomorphism 
invariance (i.e., symmetry for non-constant N^). Con- 
sequently Hfj, vanishes on shell; hence the value of the 
Hamiltonian — the conserved quantity associated with a 
local displacement N and a spatial region E — is deter- 
mined by a 2-surface integral over the region's boundary: 



E(N, E) 



H{N) 



B{N). 



(1) 



,7E 



For any choice of N this expression defines a conserved 
quasi-local quantity. Different choices of boundary term 
correspond to different boundary conditions. 

Einstein's gravity theory, general relativity (GR), can 
be formulated in several ways. For our purposes the 
most convenient is to take the orthonormal coframe t9 M = 
{>V J kdx k and the connection one-form uj a p = T a pkdx k as 
the geometric potentials. Moreover we take the connec- 
tion to be a priori metric compatible: Dg a p :— dg a p — 



0. Restricted to orthonormal frames 



where the metric components are constant, this condition 
reduces to the algebraic constraint w a ^ = uv- a ^. 

We consider the vacuum (source free) case for simplic- 
ity. GR can be obtained from the first order Lagrangian 
4-form £ GR = Q.<*P A Pa p+D^ A Tfl -V a P Aip^-^p), 
where fl a /3 '■= dio a ' p + u a 1 A uf 1 p is the curvature 2- 
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form, D-d^ :— dfi^ + uj^ v A i?" is the torsion 2-form, and 
T] aP - := *(i) a A ^ A ■ • •) is the dual form basis. The 
2-forms f2 Q/3 , U a/3 and cua are antisymmetric. We take 
k := 8ttG/c 4 — 8ir. In a "preferred boundary term" 
for GR was identified: 

B(N) = -L (A^«^ A W( / + DpN a A V J) , (2) 
107T v ' 

where A indicates the difference between the dynamic 
and reference values and Dp is the reference covariant 
derivative. The reference values can be determined by 
pullback from an embedding of the boundary into a suit- 
able reference space. 

The Energy- Extremization Program. Here we explicitly 
formulate the extremization program for static spheri- 
cally symmetric spacetimes. The Schwarzschild-like met- 
ric in "standard" spherical coordinates is given by ds 2 — 
-Adt 2 + A^dr 2 + r 2 dVtl, where dQ% = d9 2 + sin 2 9dip 2 . 
However, there are other favorable coordinate choices 
for the Schwarzschild metric (e.g., Painleve-Gullstrand, 
Eddington-Finkelstein, Kruskal-Szekeres). In order to 
accommodate most well-known coordinates, we consider 
a more general version of the Schwarzschild metric via 
a coordinate transformation t = t(u,v), r — r(u,v); the 



metric becomes ds 2 
At u t v )dudv + (A 
spacetime ds 2 — —dT 



(At 2 u -A- 1 r 2 u )du 2 + 2(A- 1 r u r v - 



1 r 2 - At 2 v )dv 2 + r 2 dQl. The Minkowski 
- - ^ 2 ' dR 2 + R 2 d0 2 + R 2 sin 2 6d$ 2 
is a natural choice for the reference. However, the es- 
sential issue of the reference choice is the identifica- 
tion between the reference and physical spacetime co- 
ordinates. A legitimate approach is to assume T = 
T(u,v),R — R(u,v),Q = 8,$> = <p and isometrically 
embed the two-sphere boundary S at (ioi r o) and its 
neighborhood into the Minkowski reference such that 
i?o := R(to,fa) = r o- Assume that the displacement vec- 
tor N = N u d u + N v d v = N l dt + N r d r = N T d T + N R d R 
is future timelike and the orientation is preserved un- 
der diffeomorphisms, i.e., \J—a :— t u r v — t v r u > and 
X- 1 := T U R V - T V R U > 0. The second term of Eq. © 
vanishes for spherically symmetric spacetimes; the energy 
can then be evaluated: 



E 

B 
C 



XT U + g 
XT, + q< 



N V C) V 3 ^, (3) 

t (R u -2r u )+g vv (R v -2r v ), (4) 
l (2r u -R u )+g uv (2r v -R v ), (5) 



where the subscripts indicate partial differentiations. 
Note that the quasi-local energy is evaluated on the 
boundary two-sphere S; the variables appearing in 
Eq. ((3|) and in the following are also evaluated on S. 
Each choice of the embedding variables {T u , T v , R u , R v } 
means a different embedding, hence a different reference. 
For any given displacement vector we extremize the en- 
ergy with respect to the embedding variables; we get four 
equations, but only three are independent: 

N U R U + N V R V = N R = 0, (6) 



X 2 T V (N U T U + N V T V ) - oT x {g uv N u + g vv N v ) = 0, (7) 
X 2 T U (N U T U + N V T V ) ~ a-^JV + g uv N v ) = 0. (8) 

A useful combination ([8|) x R v — (|7|) x R u gives 

X(N U T U + N V T V ) + a^Hg^N 11 + g vv N v )R u 

-(9uuN u + g uv N v )R v } = 0. (9) 

From Eq. © we get R u = -^R„ and N T := N U T U + 
N V T V = then R v can be found from Eq. ©: 



N u , R v , „ a(N u ) 2 

a' 1 — g(N,N) = => Rl = -4 -r- 

Rv N uyy ' ' v g(N, N) 



(10) 



We require the displacement vector to be future timelike, 
i.e., N T > and N u > 0, and the orientation to be 
preserved under diffeomorphisms, i.e., the Jacobians are 
positive. Then R v should be positive, and therefore 



R v 



a 



-N v 



Ru 



-N v 



(11) 



g(N,Ny ' — y g(N,N) 
Now we calculate the energy. Using Eq. (111[) we get 

V^a(N u B + N V C) = 2 {^/-g{N, N) - AN*\ , (12) 

where the explicit metric is used in the calculation. 
Choose N to be unit timclikc on the two-sphere, i.e., 
-1 = g(N, N) = g uu (N u ) 2 + 2g uv N u N"+g vv (N v ) 2 , then 
the quasi-local energy for any given future timelike dis- 
placement vector N reduces to 



E = r ( 1 - AN 



(13) 



which is independent of the coordinate system. The en- 
ergy expression was obtained without knowing the ex- 
plicit expressions for the variables T U ,T V . Indeed, we can- 
not solve for all four embedding variables, since there are 
only three independent equations from the extremization. 
Thus, this optimal program produces a unique energy 
with an equivalent class of references for any given phys- 
ical observer. However, it is also reasonable to impose the 
normalization condition of the displacement vector in the 
reference spacetime, i.e., g(N,N) = — 1. Using this con- 
dition and Eqs. (fflI5| fTT]> we find T u = AtuN* — A~ 1 r u N r 
and T v = At v N l - A~ x r v N r . Then the reference is 
uniquely determined. 

Now we can vary the energy with respect to the dis- 
placement vector to determine the observer who mea- 
sures the extreme energy. In view of the constraint 
g(N, N) = -A(A^ t ) 2 + yl- 1 (A r ) 2 = -1 we take y/AN 1 = 
cosh z and N r /^/A = sinh z, and then the energy extrem- 
ization dE/dz = implies coshz = 1 and QjM- \ 9e =q < 0. 
This result means that among all physical observers the 
static observer, i.e., N = A~ 1 / 2 dt, would measure the 
maximum quasi-local energy 



E = r I 



(i-VXj. 



(14) 
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For Schwarzschild this is a standard value which has been 
obtained by many researchers, e.g., 0,0,11, 12, 13j |. 

The energy expression (fTl)) can also be obtained via a 
one step approach of extremizing the energy with respect 
to the embedding variables and the displacement vector. 
We firstly require the displacement vector to be both the 
unit timelike Killing vector of the Minkowski spacctimc 



N = d T 



N u = XRa, N v = -XRu, (15) 



and a unit timelike vector in the physical spacetime, 
-1 = g(N, N), implying 

X 2 = —{g uu R v — 2g uv R u R v + g vv R 2 ) , (16) 

With these assumptions, the energy © reduces to 



E = -X 
2 



1-q 1 X 2 + 2a 1 (g uv r u R v - g uu r v R v 



(17) 



Now there are only two embedding variables appearing 
in the energy expression. Varying the above energy ex- 
pression with respect to these two variables we get 

{gwRu - g U vRv){i + a~ l X~ 2 ) 

+2R v (r u R v -r v R u ) = 0, (18) 

(guuRv - sw-R«)(l + or x X~ 2 ) 

-2R u (r u R v -r v R u ) = 0. (19) 

The combination (TT8")) x i?„ + p^|) x R v gives X^ 1 = y/—a; 
then both (jTSj) and l[T9)) reduce to the condition r u R v — 
r v R u = 0. Together with Eq. ^ we get R* = r 2 u /A 
and R 2 — r 2 /A. We should pick the plus sign for R u and 
R v . The unique energy produced by this program is then 
E = r(l — V^4), which agrees with what we found above. 

Energy Measured by Various Observers. For the 
Schwarzschild spacctimc A = l — 2m/r; it is obvious that 
Eq. (fl4|) can only be valid outside the black hole horizon: 
there is no static observer inside the black hole. In order 
to discuss the energy inside a black hole, let us examine 
our energy formula Eq. (|13j) for a radial geodesic observer 
in the Schwarzschild spacetime. For an observer who falls 
initially with velocity vq from a constant distance r = a, 
there are two different types of orbits: (1) the crash 
orbit — where ingoing observers crash directly into the 
singularity at r — 0, while outgoing observers first shoot 
out to the turning point r max and then fall back and crash 
into the singularity; (2) the crash/escape orbit — where 
ingoing observers crash, but outgoing observers can es- 
cape to infinity by having a large enough initial velocity. 
The displacement vector for such observers is the unit 
tangent of the geodesic; then N l = 1 _ 2 1 m /r " ' 

where 2m < a and < r < r max . The energy mea- 
sured by this observer is E = r (l — ^/ ^-T^" ") ■ 



result agrees with that of the Brown- York quasi- local en- 
ergy expression 1J, |15J. One can see that the energy 
decreases as the initial velocity increases. When the 
initial velocity Vo is less, equal, or greater than ^J2mja 
(which is the escape velocity from the Newtonian point 
of view) the energy is positive, zero, or negative, respec- 
tively. The negative value for the energy may appear 
odd, but it can be explained physically. It is correlated 
with the geometric property that the scalar curvature 
of the spacelike hypersurface orthogonal to the displace- 
ment vector is (unlike the usual cases) negative. Note 
that the ingoing geodesic observers can measure energy 
inside the black hole. This energy is proportional to the 
radial distance r, so it is a smooth function in the region 
< r < r max . 

In addition to the static and radial geodesic ob- 
servers there are other natural choices — in particular 
the unit normal of the constant coordinate time hy- 
persurface in various coordinate systems. Using this 
idea in ingoing Eddington-Finkelstein coordinates, ds 2 = 

r 2 dn?,, where 
r[l - 



2/ndu 



T-\J l-\-2m/r 

3pA _1 dr, dv 

-1/21 



l+2m/r ' 

du = dt + 

(1 + 2m/r)~ 
Whereas for 
ds 2 — —du 2 

A^ 1 ^J2m/rdr, dv 



+ y/l + 2m/rdv 



dr, gives -Eef 
2m/r 



dt 



2m(l + 2m/r + yjl + 2m /r) . 
ingoing Painleve-Gullstrand coordinates, 
( y ^2m/rdu + dv) + r 2 dfl2: where du = 
dr, the energy vanishes 



All 



these outcomes are smoothly dependent on r. 



The Geometrical Meaning. The proposed optimal pro- 
gram can be interpreted as an adapted coordinate choice 
for any given displacement vector, i.e., observer. The 
adapted coordinates, denoted as {u',v'}, have the asso- 
ciated coframes 

i9° = aP u ,du', d 1 = a\,dv! + a^dv 1 , (20) 

and the displacement vector is the unit timelike vector, 
namely N = eo; its components can be expressed in 
terms of (a°, : a*,, a\,). The embedding variables are de- 
termined by requiring on S 

i9° = a a a ,du' = dT = T u ,du' + T v ,dv', 
i? 1 = a\ L ,du! + a l u ,dv' = dR = R u ,du! + R v ,dv'. (21) 



This 



This identification leads directly to (| 13|) . Thus, the 
extremization program, in these adapted coordinates, 
yields an isometric embedding on the two-sphere bound- 
ary in a way such that N is the unit timelike vector in 
both the physical and Minkowski reference spacetimes. 

Dynamic Spherically Symmetric Spacetimes. The opti- 
mal reference choice program can likewise be applied to 
dynamic spherically symmetric spacetimes, the FLRW 
cosmological models. With a calculation similar to that 
above, we found that the optimal energy for a given dis- 
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placement vector is 



E — ar 



1 - kr 2 N l 



VI - kr 2 



N r ), (22) 



where a = da/dt and k = —1, 0, 1 is the sign of the spatial 
curvature. An obvious choice is the comoving observer, 
N = dt, then the energy is E = ar (l — yl — fcr 2 ) — a 
value which has been found previously [16| — which is neg- 
ative for negative curvature. On the other hand, varying 
the energy ((22)) with respect to the displacement vector 
we find the maximum energy 



E = 



ar 3 (k + a 2 ) 



1 + \/T— kr 2 — a 2 r 2 ' 



(23) 



(which again does not exist for all observer locations r). 
The associated displacement vector is 



TV 



Vl - kr 2 



\Jl — kr 2 — a 2 r 2 



ar Vl - kr 2 

dt n i 2 -2 =7^' ( 24 ) 

a yl - fcr z — a z r z 



which is just the unit dual mean curvature vector 17 1. 
For a dust model with matter energy density p, by im- 
posing the Friedmann equation, k + a 2 — (8ir/3)pa 2 , the 
energy expression (f2"3"| becomes 



E = 



(8Tr/3)pa 3 r 3 



1 + y/l - (87r/3)pa 2 r 2 



> 0. 



(25) 



For small r the proper matter interior limit required by 
the equivalence principle is evident. While we find that 
the observer dependent quasi-local energy can be nega- 
tive, for the particular observer who measures the maxi- 
mum value, the quasi-local energy is positive. 

Conclusion. The covariant Hamiltonian quasi-local en- 
ergy expression @ has certain nice properties, @, H, 0, 

El 



111 ], but it suffers from two ambiguities: which dis- 
placement vector and which reference? We have pro- 
posed isometrically embedding the two-sphere boundary 
and its neighborhood in the dynamic spacetime into the 
Minkowski reference and then extremizing the energy 
to determine the embedding variables. Here we have 
discussed the application of this program to static and 
dynamic spherically symmetric spacetimes where we ob- 
tain, for any given future timelike displacement, a unique 
quasi-local energy which can be positive, zero, or even 
negative. When we further vary the energy with re- 
spect to the displacement vector we can identify a spe- 
cial observer who measures the maximum — and always 
positive — energy. 

Moreover, the optimal program is, at least in this 
spherical case, actually equivalent to an adapted coor- 
dinate choice for each observer. In such coordinates, 



our program not only isometrically matches the geome- 
try near the two-sphere boundary, but also identifies the 
displacement vector as the unit timelike vector orthog- 
onal to the constant coordinate time hypersurface, and, 
at the same time, as the unit timelike Killing vector of 
the Minkowski reference. This observation lends further 
support to the proposed optimal program. 

We believe that this spherical case is the main test case; 
it shows that our program has promise as a universal 
approach for determining the reference needed for the 
covariant Hamiltonian boundary term quasi-local energy 
for general spacetimes. 
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